Quasi-Normal Modes of a Schwarzschild White Hole 



Nigel T. Bishop and Amos S. Kubeka 
Department of Mathematical Sciences, University of South Africa, P. O. Box 392, Unisa 0003, South Africa 

We investigate perturbations of the Schwarzschild geometry using a linearization of the Einstein 
vacuum equations within a Bondi-Sachs, or null cone, formalism. We develop a numerical method 
to calculate the quasi-normal modes, and present results for the case £ — 2. The values obtained 
are different to those of a Schwarzschild black hole, and we interpret them as quasi-normal modes 
of a Schwarzschild white hole. 



I. INTRODUCTION 
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The theory of a linear perturbation of a black hole was 
developed some time ago 111, 0, 0, 0); see also the text- 
book p] and the review p|. The essential idea is that 
the vacuum Einstein equations are linearized about the 
Schwarzschild (or Kerr) geometry described by the usual 
(i, r, 8, <fi) coordinates. Then a standard separation of 
variables ansatz is applied, with metric quantities behav- 
ing as an unknown function of r y,Yp m {9, <fi) exp(ii>t) (ac- 
tually, the angular dependence is somewhat more com- 
plicated, and the technical details can be found in the 
literature). There results an ordinary differential equa- 
tion in r, and the quasi-normal modes are obtained by 
finding the special values of v for which solutions ex- 
ist that satisfy appropriate boundary conditions in the 
neighbouhood of the event horizon, and of infinity. 

Quasi-normal mode theory has become a cornerstone 
of modern general relativity theory. They have been seen 
in numerical relativity simulations of binary black hole 
coalescence. And, while not yet actually observed, it 
is strongly expected that they will be measured by the 
LIGO collaboration, and certainly by LISA, yielding pre- 
cise information about the parameters describing a black 
hole from some coalescence event. 

In the usual approach to linear perturbations of a 
black hole, the linearization is performed using standard 
Schwarzschild (or Kerr) coordinates (t,r,6,(/>). It is also 
possible to perform the linearization using Bondi-Sachs 
coordinates, which is a coordinate system based on out- 
going null cones. This has been done in previous work, in 
order to obtain analytic solutions of the linearized Ein- 
stein equations for the purpose of testing numerical rel- 
ativity codes. As with the usual approach, one ends up 
with a second order ordinary differential equation involv- 
ing I and v as parameters, Eq. ([5]). However, when the 
quasi-normal modes were calculated for this equation, it 
was found that they are not the standard ones. Differ- 
ent physical problems are considered in the two cases, 
as illustrated in the Penrose diagram of Schwarzschild 
spacetime (Fig. [I}. K is a typical hypersurface used in 
finding the quasi-normal modes of a black hole, and the 
direction of wave propagation at the boundaries of K is 
shown by arrows. On the other hand, N is a typical 
hypersurface used in finding the quasi-normal modes of 
Eq. (0. From the direction of wave propagation on N, 
the resulting quasi-normal modes can be interpreted as 




FIG. 1: Penrose diagram illustrating the differences, in terms 
of location and boundary conditions, between the hypersur- 
faces K and N. 

The plan of this paper is as follows. Section [II] sum- 
marizes previous work on the Bondi-Sachs metric and 
linearized solutions within that framework. Section IIIII 
describes our approach to calculating the quasi-normal 
modes, and Sec. IIVI presents the results. We end with a 
Conclusion, Sec. IVl 



II. BACKGROUND MATERIAL 

The Bondi-Sachs formalism uses coordinates x l = 
(u, r, x A ) based upon a family of outgoing null hyper- 
surfaces. We label the hypersurfaces by u =constant, 
null rays by x A (A = 2,3), and the surface area coor- 
dinate by r. In this coordinate system the Bondi-Sachs 
metric 0, 0, 0] takes the form 



ds 1 = - 



,2/J 



l+W)-r 2 h AB U A U B 



du 2 - 2e 2l3 dudr 



2r z h AB U"dudx A + r 1 h ab^ dx 



(1) 



where h AB h,Bc = $b an< ^ det(fiAB) = det{q A B), with 
qAB being a unit sphere metric. We represent qAB by 
means of a complex dyad qa- For example, in the case 
that the angular coordinates are spherical polar (0,<f>), 
the dyad takes the form 



q A = (l,isin0). 



(2) 
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For an arbitrary Bondi-Sachs metric, Hab can be repre- 
sented by its dyad component 



J = h AB q A q B /2, 



(3) 



We also introduce the spin- weighted field U = U A qA, as 
well as the (complex differential) eth operators 3 and 3 
(Tot - In Schwarzschild space-time, W = —2M, (3 = 0, 
U A = and J = 

We use Zi m , rather than Yg, m , as spherical harmonic 
basis functions, where [Hj | the have orthonormal 
properties similar to those of the Y( m , and are real. 

We assume the following ansatz, representing a small 
perturbation of the Schwarzschild geometry 



(3 = ^{Mry uu )Z lm , U = ^{U {r)e l 
J = Vt(J (r)e iuu )<5 2 Z em , 
W = -2M + R(w (r)e ivv )Z em . 



(4) 



Using the above ansatz, Ref. [ll| constructed the re- 
sulting linearized Einstein vacuum equations. As ex- 
pected, the angular and time dependence factored out, 
and a system of ordinary differential equations (in r) was 
obtained. As discussed in Ref. the system can be 
manipulated to give 



x 3 (l - 2xM)^- + 2^-(2x 2 + ivx - 7x 3 M) 
ax ax 

-2(x(f +1- 2)/2 + 8Mx 2 + iv)J 2 = 



(5) 



where Ji{x) = d 2 J /dx 2 and x = 1/r. (Actually, 
Ref. [ll| gave Eq. © only in the case 1 = 2, and here 
we give the formula for general I). 



III. PROBLEM SPECIFICATION 

We note that Eq. ([5|) has singularities at x = and 
x = 0.5 M. The problem is to find values of v for which 
there exists a solution to Eq. |(SJ) that is regular every- 
where in the interval [0, 0.5M]; these values of v are the 
quasi-normal modes. This is the same situation that is 
faced when finding the quasi-normal modes of a black 
hole. The first solution to this problem was obtained 
by using series solutions around the singular points, and 
a numerical solution of an ordinary differential equation 
within the interior of the interval Subsequently, it 
was shown [l2[ how the theory of 3-term recurrence re- 
lations [lH for the series solution about x = 0.5Af could 
be used to determine the quasi-normal modes. 

It is straightforward to write Eq. ([5]) with the origin 
transferred to x — 0.5M , and then to evaluate the recur- 
rence relation satisfied by a regular solution (see Eqs. (j2"4"|) 
and (|2"5|) below) . We find a 4-term, rather than a 3-term, 
recurrence relation. While it may be that the quasi- 
normal modes could be found by an approach similar to 
that of [HJ , this is not a practical option since there does 
not seem to be available a well-developed mathematical 
theory of 4-term recurrence relations. 



Instead, we proceed along the lines used in We 
construct the asymptotic series about the essential sin- 
gularity at x — 0, and use it to find a solution to within 
a specified tolerance at a point xo > 0. We then use this 
solution as initial data for a numerical solution of Eq. J5|) 
in the range (xo, x c ) where x c < 0.5M; actually, as in [4|, 
we do not integrate Eq. directly but first convert it 
to first-order Ricatti form. Finally, we construct the reg- 
ular series solution about x = 0.5M and use it to find a 
solution at x — x c . Then a value of v is a quasi-normal 
mode if the difference at x = x c between the regular 
series solution and the numerical solution, vanishes. 



A. Asymptotic series solution about the essential 
singularity at x = 

Since the singularity is essential, the resulting scries 
solution has radius of convergence zero, although it is 
asymptotic. We use [l4| to determine rigorous bounds 
on the error of approximating the solution by its first n 
terms. Note that a series solution Jzix) = XmLi a nX n to 
Eq. ([5]) can be generated by the recurrence relation 

n 2 + n — 6 2n(n + 2) 
a n = -a n -i-zT—, — + a ra _ 2 M — , (6) 



2iv{n - 1) 



2iv(n - 1) 



with ai = 1, a,2 = 0. 

In order to use the theory developed in [lij], we must 
first transform Eq. ([5]) to asymptotic form by 



1 



x — > z = — , 

x 



(7) 



and investigate the solution about the singularity at in- 
finity. We find 

z 2 (z - 2) — f^- - z(2z 2 iv + 2z- 10) — 
dz 2 dz 

-{2z 2 w + ^z+ 16)J 2 (z) = 0, (8) 

where we have normalized the scaling of z by setting 
M = 1. We evaluate quantities used in [Til ]: 

2{z-b + ivz 2 ) _ 2(2z + 8 + ivz 2 ) 
* ~ z(z - 2) ' 9 ~ z 2 (z- 2) ' 

f = -2iv, f x = -2- Uu, g Q = 0, 
gi = —2iv, p = if, cr = 2 + 2iv. (9) 

Then the solutions can be written as 



J 2j (z)^exp(\ j z)z^J2^T 

s=0 



(10) 



where 

Ai = 0, mi = -1, A 2 = 2iu, ^2 = 3 + tiu. (11) 
Following [l4]], we let the solution to Eq. © be 

Ja(*) =L n (z) + e n (z) (12) 
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where L n (z) = exp(\\z)z^ 1 — — 

Z — l z s 



(13) 



and define the residual R n (z) by 
d 2 L n (z) sdL n {z) 



dz 2 



•/(* 



dz 



+ g(z)L n (z) = 



R n {z) 



with 



(14) 



(15) 



in some region \z\ > 6 and where B n is calculable. 
Ref. [lij obtains a bound on e n (z) provided the quan- 
tity C(n, 6, v) defined immediately below satisfies C < 1, 
where 



C(n, 6, v) 



{3^Y(\(n + !) + !)) 
|2w|r(|(n + l) + i)) (n + 1) 

where /3 is bounded by 



(16) 



/3 < |4ii/| + 



6-2 



32 



1 



6(6-2) 





3-4iv 




f|2 + 4»V| + 




) 




6-2 





(17) 



Given v and 6, we use numerics to determine conditions 
on n such that C < 0.99 and then we bound e n (z) by 



2£L 



(18) 



/?(1-C(n,6, 

We also need to bound the error e^(x) in using a finite 



series to estimate dJ j[ x> . Noting that 

2 dJ 2 (z) 



dx 

dJ 2 (x) 
dx 



= — z 



dz 



the bound on the error is 

K(x)| < 



2\iv\B„ 



P(l-C(n,b, u)\z\ 



n-l ' 



(19) 



(20) 



1. Numerical implementation 

We have written Matlab code that takes as input v 
and 6, and then finds /3 and the lowest value of n such 
that C < 0.99. Then the code finds the maximum of 
the absolute values of e n (b) and e' n (x = 1/6). A bisec- 
tion method program takes v as input and refines 6 until 
the absolute value of the maximum error is in the range 
(0.5, l)x machine precision (about 2 x 10 -16 ). The code 
returns the values of 1/6 and L n (b) / L' n (x = 1/6). 



B. Numerical integration of Eq. ([5]) 

The first step is to transform Eq. ([5]) into first-order 
Ricatti form. Defining a new dependent variable v(x) by 



we obtain 



^2(x) — > v(x) 



dv 



1 dJ 2 (x) 
J 2 (x) dx 



(21) 



x 3 (l - 2x) — +v 2 
\dx 

+ 2x(2x +iv- 7x 2 )v - 2(2ir + 8x 2 + iv) = 0. (22) 



The numerical integration of Eq. (|22p near the singu- 
larity at x = can be tricky, because we need the result 
to be as accurate as possible. We found that a fourth 
order Runge-Kutta scheme (ode45 in Matlab) performed 
better than the stiff schemes, provided stringent toler- 
ance conditions were used (specifically, RelTol = 10 -12 , 
AbsTol = 10~ 12 , MaxStep = 2 x 10" 6 ). Under these 
conditions each integration to x c (=0.25) takes of order 
100s. 



C. Series solution about the regular singularity at 

x = 0.5M 



We first make the transformation 
x — > s — 1 — 2x 
to Eq. (O and obtain 



s (i- s y 



d 2 J 2 (s) 



ds 2 

- 4(m/ + 3 - 5s + 2s 2 ) J 2 (s) = 0. 



-(l-s)(4w-3 + 10s-7s 2 ) 



(23) 



dJ 2 (s) 
ds 
(24) 



This equation has a series solution X)o° a nS n that satisfies 
the recurrence relation 



3 + iv 

ao = 1, ffli = 4- — , a 2 



15(4 + 3iv) 



3-Aiv J * 2(1 - iv)(3 - Aiv) 

4niv — 8iv — 5 — in 2 — An 

1 n(Aiv — n — 2) 

4 + 3n 2 + 2n (l-n)(l + n) 
+ a«-2— - + a»- 3 lA . ^"-(25) 



' n(Aiv — n — 2) 



' n(Aiv — n — 2) 



The radius of convergence of the above series is s < 1, 
and, given ^, the numerical evaluation of the coefficients, 
and then of the series, is straightforward. Using x c = 
0.25 means that we need to evaluate the series at s = 
0.5. We terminate summation of the series at the first 
term smaller than 10 -18 (typically, about 60 terms), and 
thus expect the result to be accurate to within machine 
precision (about 2 x 10~ 16 ). 
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IV. RESULTS 

We have written a Matlab program that, given a value 
of v, first uses the asymptotic series to find the value 
v of v(x) (as defined in Eq. (|2Tjl ) at x = x = 1/b, 
and then integrates numerically Eq. ([22]) between Xq and 
x c = 0.25, obtaining a complex number v + = v(x c ); and 
secondly uses the regular series about x — 0.5 to find 
V- = v(x c ). Defining 

9u=v+-v-, (26) 

the quasi-normal modes are those values of v such that 
g v is indistinguishable from zero. 

We calculated g v for values of v in the range v — a + ib, 
0.1 < a < 1.07, 0.05 < b < 0.89, in increments of 0.03. 
The results are shown in the contour plot in Fig [21 The 
green line is the zero contour of 3(g„), the red line is the 
zero contour of 3fJ(<?„), and the blue line is the boundary 
of a region where the computation is probably unreliable 
(because the computed curve oscillates, indicating that 
a smaller step- length is required). Clearly, the quasi- 
normal modes lie at the intersection of a red and a green 
line, and from the plot we can read off an estimate for 
the lowest mode, v = 0.9 + 0.63i. We then applied a 
secant method, obtaining a final estimate for the lowest 
quasi-normal mode at 

^ = 0.883 + 0.614i. (27) 

In this case, x = 0.036493228795438, v = 
0.036838521818950 + 0.000637428772012i, (3 
0.988517790240599, and 62 terms were used in the 
asymptotic series. The contour plot indicates another 
quasi-normal mode at about v = 1.06 + 0.63z, but we 
did not investigate further. 




FIG. 2: Contour plot in the complex plane of v showing the 
contours where 9f(<7„) = (red) and Ss(g v ) = (green). 

We now use the value in Eq. (|27[) . and vary the nu- 
merical methods so as to determine the accuracy with 



which g v has been determined. In Fig [3] the integration 
between xq and x c is carried out with different values of 
MaxStep, 2 x 10~ 6 , 10~ 6 and 5 x 10~ 7 , and also an error 
of an amount (1+z) x 10 -15 is introduced into the value of 
vo at xq in the case MaxStep = 2 x 10~ 6 . Also, numerical 
integration of Eq. ([2"2")l as well as a series solution is used 
in the range (x c ,0.5). The various curves lie on top of 
each other and are visually indistinguishable. Taking all 
these options into account, the maximum value noted for 
g v was (6.02 + 5.87i) x 10 -4 . Using intermediate results 
from the secant root-finding process to estimate 

— = 3.95 + 0.69i, (28) 

it follows that the possible error in Eq. ([2T| is 

|(3.95 + 0.69£) x (6.02 + 5.87?) x 10 4 | = 0.003, (29) 
so that Eq. ([2~7| should be amended to read 

v = 0.883 + 0.614i + 0.003A; (30) 
where k is a complex number satisfying \k\ < 1. 

0.3 1 1 1 1 1 1 1 1 1 1 1 
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FIG. 3: The real (solid line) and imaginary (dotted line) parts 
of v(x) in the quasi-normal mode case v = 0.883 + 0.614i. 

The lowest quasi- normal mode of a Schwarzschild black 
hole is at v — 0. 37367+0. 08896i. We have used this value 
in our code, and obtained Fig [H from which it is clear 
that this value of v is not a quasi- normal mode of Eq. ([5]). 

V. CONCLUSION 

Using a linearization of the vacuum Einstein equations 
about the Schwarzschild geometry, within a Bondi-Sachs 
framework, we have constructed a numerical procedure 
to calculate the quasi-normal modes. The value of the 
lowest mode in the case t = 2 is not a quasi-normal 
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FIG. 4: The real (solid line) and imaginary (dotted line) parts 
of v(x) in the case v = 0. 37367 + 0. 08896i, indicating that the 
lowest quasi-normal mode of a Schwarzschild black hole is not 
a quasi-normal mode of Eq. ([5]). 



Eq. ([5]), and thus it is important to discuss the extent to 
which this has been verified. Eq. ([5]) was derived in [ill ], 
and there has been no subsequent, independent, deriva- 
tion. Nevertheless, Eq. ([5]) has been subject to some con- 
sistency checks, since ref. confirmed that solutions 
obtained also satisfy the remaining Einstein equations 
(the constraint equations). Further, in the case M = 0, 
solutions based on Eq. |5]| have been used as analytic so- 
lutions for the testing of numerical relativity codes based 
on the Bondi-Sachs metric, and the expected order of 
convergence was observed [H, [l6| . 

The evidence for the existence of black holes is now 
very strong. However, the question about the existence 
of white holes is much more problematic, since such ob- 
jects cannot form from regular initial data, but instead 
must have been created as part of the creation of the 
universe. The present work provides a possible observa- 
tional signature of a white hole, since it is in principle 
possible for a gravitational wave detector to extract the 
parameters of a quasi-normal mode from a gravitational 
wave signal. 



mode of a Schwarzschild black hole, and further the low- 
est quasi-normal mode of a Schwarzschild black hole is 
not a quasi-normal mode of Eq. As discussed in the 
Introduction, this apparent discrepancy can be avoided 
by interpreting the quasi-normal modes of Eq. ([5]) as be- 
ing those of a white hole rather than those of a black 
hole. 

The results obtained depend crucially on the validity of 
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